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Abstract 

In this paper, we propose & implement on WMAP 5-year data, a model independent approach 
of foreground power spectrum estimation for multifrequency observations of CMB experiments. 
Recently a model independent approach of CMB power spectrum estimation was proposed by Saha 
et al. 2006. This methodology demonstrates that CMB power spectrum can be reliably estimated 
solely from WMAP data without assuming any template models for the foreground components. In 
the current paper, we extend this work to estimate the galactic foreground power spectrum using the 
WMAP 5 year maps following a self contained analysis. We apply the model independent method in 
harmonic basis to estimate the foreground power spectrum and frequency dependence of combined 
foregrounds. We also study the behaviour of synchrotron spectral index variation over different 
regions of the sky. We compare our results with those obtained from MEM foreground maps which 
are formed in pixel space. We find that relative to our model independent estimates MEM maps 
overestimates the foreground power close to galactic plane and underestimates it at high latitudes. 



I. INTRODUCTION 

The study of Galactic foreground emission in the range from few MHz to few hundreds of GHz is very important 
for CMB observations. Characterising the foreground emission unravels the galactic physics. Its understanding leads 
to more reliable removal of foreground contamination from the CMB temperature or polarization anisotropy maps. 
Diffuse galactic foregrounds consist of three main components : dust, free- free and synchrotron emission. Recent 
literature also report evidence of other possible contaminating components like spinning dust, hard synchrotron and 
galactic haze at WMAP frequencies [5, 6, 16, 17]. It is important to measure the galactic foregrounds components 
and their spectral behaviour unbiased by prior expectations. The physical mechanism and spectral behaviour of three 
main foreground components are summarised below: 

• Synchrotron emission arises when an electron moves at a relativistic velocity along a magnetic field line. In 
terms of antenna temperature the frequency dependence of synchrotron emission can be written as, 

T s (x iA , (1) 

where the spectral index, /3 S , varies across the sky. Based on the WMAP data, WMAP team estimated /3 S ~ —3.5 
at high latitudes and f3 s ~ —2.5 close to star forming regions near the galactic plane. 

• Free-free emission arises when free electrons passing through the hot interstellar medium are deflected and 
slowed down by ionized atoms mostly protons. The change of kinetic energy due to deflection is converted into 
free-free emission. The frequency dependence in terms of antenna temperature can be written as, 

Tf oc zA , (2) 

where the index /3f is flatter compared to f3 s and can be approximated as /3f = —2.14 in the WMAP frequency 
range. 

• Dust emission occurs when dust grains in the interstellar medium heated by the photon flux from the stars seek 
thermal equilibrium by emitting in infra-red or far infra-red range of the spectrum. Dust emission can be well 
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modeled by grey body spectrum and the frequency dependence is given by, 



exp{hv / KT dust ) - 1 



(his « KT dust ) . 



(3) 



Different techniques, such as, Maximum Entropy Method (MEM) [4, 6, 8], Correlated Component Analysis (CCA) 
[10], have been studied in the literature for modelling individual foreground components in pixel space. However, 
all these methods require foreground models in terms of templates, e.g., 408MHz all sky synchrotron emission map 
(Haslam et al. 1982) [13], dust map at 94GHz by Finkbeiner, Davis & Schlegel (FDS99) [15] and Finkbeiner Halpha 
maps [16]. While modeling foreground components these methods rely on extrapolation of the external templates to 
WMAP frequencies. There is always an uncertainty in extrapolating low frequency or high frequency template to 
WMAP frequencies [11]. In this paper, we discuss about a model independent estimation of foreground behaviour 
without using any information from external template. 

A model independent foreground removal method from the multifrequency CMB data was first proposed by Tegmark 
& Efstathiou 1996, and was implemented on the WMAP 1-year data by Tegmark et al. 2003. This method was 
extended by [1, 2, 3] to estimate CMB power spectrum by cross-correlating cleaned maps with independent detector 
noise. We extend this method to get the spatial and spectral distributions of diffuse foreground power spectrum. 
The main advantage of this method is that it does not require any assumptions about the foreground emission. It 
is based only on the fact that CMB anisotropy is independent of frequency in thermodynamic temperature unit 
whereas foregrounds have a frequency dependence. In this paper, we motivate ourselves to estimate power spectrum 
of the composite emission due to all diffuse foreground components in a self-contained manner using WMAP data 
only. However, we do not address the problem to present power spectrum due to individual foreground components. 
The method does not require extrapolation of foreground templates from outside frequency or fitting any foreground 
template. Amongst other advantages are, it is computationally fast and unlike MEM or CCA our composite foreground 
maps are not limited by the lowest resolution frequency band. By extending the method we provide a map of the 
variation of the synchrotron spectral index over different regions of the sky. 

The rest of the paper is organized as follows. Section II briefly describes the model independent CMB power 
spectrum estimation methodology. Section III describes in detail the model independent foreground power spectrum 
estimation method. Section IV is dedicated to the estimation of synchrotron spectral index over different regions of 
the sky. Finally we conclude in Section V. 



The basic assumption of model independent CMB power spectrum estimation is that CMB contributes equally in 
all frequency channels (in terms of thermodynamic temperature). The temperature anisotropy at a frequency channel 
i can be written as, 



where AT C (n) and ATf(n) are respectively the CMB and foreground components and AT- Ar (n) is the noise compo- 
nent. Here Bi(n.n') is the circularly symmetric beam denoting the smoothing of the map due to the finite resolution 
of the detector. The spherical harmonic transform of the map at the frequency channel i can be written as, 



where the index i runs from 1 to 5 for five frequency bands of WMAP. In our analysis, we perform foreground removal 
on harmonic space and separately over different regions of the sky. We first partition the entire sky in nine different 
regions according to their level of foreground contamination as described in [2, 7] and then perform foreground removal 
for each region iteratively. Below we briefly describe the basic algorithm of the procedure. 
We define a cleaned map as a sum of linearly weighted 5 WMAP maps as follows, 



II. REVIEW OF THE CMB POWER SPECTRUM ESTIMATION 




(4) 
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(5) 



where the weight factor depends on the frequency i as well as on the multipole I. Since the frequency channels of 
WMAP are of different resolution, we deconvolve each map by the corresponding beam, B\, prior to linear combination. 
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The power spectrum of the cleaned map is given by, 

1 1 

r~iclean \ ^ „clean „ 

c i - + i ^ lm 



clean* 
lm 



21 + 1 2^ 2^ i B f i B i ■ 

The above equation can be simplified and written as a matrix equation as, 

C clean = WlCl Wf 

= W,(Cf + Cf)Wf 

= C^W L ee T Wf + W^Cf Wf, 

where e = (1, 1, ..,1) is a column vector of five ones, 'Cf ' and 'C^' represent respectively the CMB and non-CMB 
( i.e., foreground plus detector noise) covariance matrices. 'Cf ' is the CMB power spectrum. To preserve the CMB 
power spectrum in the cleaned map we impose the condition, W^e = e T Wf = 1. Hence we can re-express above 
equation as follows, 

cf ean = cf + w z cfw^. 

Since in the above equation, CMB power spectrum remains independent on weights, minimizing Ci lean minimizes 
the combined contamination coming from foregrounds and the detector noise without altering the CMB power. 
Minimization of Cf ean = W/C/W^ with the condition that W/e = e T W^ = 1 is a standard Lagrangian multiplier 
problem and has a well known solution, 



cr x e 

-i — 



Putting this back in equation (5), we get the cleaned CMB map. The cleaned map in fourier space can be written 

as, 

af! n an (i)=af m (i)+af^(i), (7) 
where 'RN' denotes the residual noise in the cleaned map. 



III. ESTIMATION OF FOREGROUND POWER SPECTRUM 



The WMAP satellite has 10 differential assemblies distributed over five different frequencies ranging from 23GHz 
to 94GHz. K and Ka band have one differential assembly (DA) each. Q and V band have two DA's namely Ql, Q2, 
VI and V2 respectively. W band has four DA's namely Wl, W2, W3 and W4. For W band, we simply average the 
pairs of DA's in the W band to form 6 DA maps - W12, W13, W14, W23, W24, W34 out of 4 DA's of W band. In 
our analysis, we smooth all the 5Yr WMAP maps to a common beam resolution of 1°. To estimate the cleaned CMB 
maps, we linearly combine the 4 DA's of different frequencies at time as described in Sec II and in details in [2, 3]. 
The various different 4 channel combination and 3 channel combination cleaned maps that can be obtained are listed 
in table I. 

The spherical harmonic transform of the WMAP maps at each frequency channel referred by index z, smoothed at 
1° beam resolution can be written as, 

a\m = B i a ?m + B ia(m(i) + Q>hn(i)- 
The one degree smooth cleaned map can be written as, 

„clean d r ,C , td ,,RN /q\ 

a lm = B l a lm + B l a lm > \P) 

where 'RN 1 denotes the residual noise in the single cleaned map obtained by our model independent analysis and B\ 
is the fourier transform of the beam at one degree resolution. To obtain the foreground power spectrum, we subtract 



4 



4-channel combinations (n c = 4) 



(K,KA)+Q1+V1+W12=(C1,CA1) 
(K,KA)+Q1+V1+W13=(C2,CA2) 
(K,KA)+Q1+V1+W14=(C3,CA3) 
(K,KA)+Q1+V1+W23=(C4,CA4) 
(K,KA)+Q1+V1+W24=(C5,CA5) 
(K,KA)+Q1+V1+W34=(C6,CA6) 
(K,KA)+Q2+V2+W12=(C7,CA7) 
(K,KA)+Q2+V2+W13=(C8,CA8) 
(K,KA)+Q2+V2+W14=(C9,CA9) 
(K,KA)+Q2+V2+W23=(C10,CA10) 
(K,KA)+Q2+V2+W24=(C11,CA11) 
(K,KA)+Q2+V2+W34=(C12,CA12) 



(K,KA)+Q1+V2+W12= 
(K,KA)+Q1+V2+W13= 
(K,KA)+Q1+V2+W14= 
(K,KA)+Q1+V2+W23= 
(K,KA)+Q1+V2+W24= 
(K,KA)+Q1+V2+W34= 
(K,KA)+Q2+V1+W12= 
(K,KA)+Q2+V1+W13= 
(K,KA)+Q2+V1+W14= 
(K,KA)+Q2+V1+W23= 
(K,KA)+Q2+V1+W24= 
(K,KA)+Q2+V1+W34= 



:(C13,CA13) 
:(C14,CA14) 
: (C15,CA15) 
: (C16,CA16) 
:(C17,CA17) 
: (C18,CA18) 
:(C19,CA19) 
: (C20,CA20) 
: (C21,CA21) 
: (C22,CA22) 
: (C23,CA23) 
: (C24,CA24) 



3-channel combinations (n c = 


3) 


Q1+V1+W12= 


=C1 


Ql+V2+W12= 


=C13 


Q1+V1+W13= 


=C2 


Ql+V2+W13= 


=C14 


Q1+V1+W14= 


=C3 


Ql+V2+W14= 


=C15 


Q1+V1+W23= 


=C4 


Ql+V2+W23= 


=C16 


Q1+V1+W24= 


=C5 


Ql+V2+W24= 


--C17 


Q1+V1+W34= 


=C6 


Ql+V2+W34= 


--C18 


Q2+V2+W12= 


=C7 


Q2+V1+W12= 


--C19 


Q2+V2+W13= 


=C8 


Q2+V1+W13= 


--C20 


Q2+V2+W14= 


=C9 


Q2+V1+W14= 


=C21 


Q2+V2+W23= 


=C10 


Q2+V1+W23= 


=C22 


Q2+V2+W24= 


--C11 


Q2+V1+W24= 


=C23 


Q2+V2+W34= 


--C12 


Q2+V1+W34= 


=C24 



TABLE I: The table shows 48 different combinations of the DA maps used in our 4 channel cleaning method and list of the 24 
possible combinations in the 3 channel cleaning method [2, 3]. 



the cleaned maps from the DA maps obtained as described earlier. The spherical harmonic transform of Di's (see 
table II) which is obtained after subtracting the cleaned cmb maps at each frequency channel can now be written as, 

4n - *?™ n = B ia f m (i) + (a? m (i) - B ia ™). 
= Biaf m (i) +a^(i). 

To estimate the foreground power spectrum at 1° beam resolution, we remove the noise bias by cross correlating pairs 
of CMB free maps which have no common DA/detector in the cleaned maps involved. Assuming there is no cross 
correlation between the foreground and noise and independent noises for two different detector, we can write the 1° 
beam smoothed foreground power spectrum Cf as, 

{{aJZ ta \i) - at an (i))(a™ al *(j) - af™ n *(j))) 

= (B*af m {i)af*{j)) + (a» m (i)aC(j)). 

= Cf. (9) 



CMB subtracted WMAP Q band 


Cross combinations 


Ql- 


C 


(K,KA)+Q1+V1+W12] 


= (D01,DA01) 


Q2- 


C 


(K,KA)+Q2+V2+W34] = 


(D13,DA13) 


D01 d 


$ DA13 


DA01 d 


$ D13 


Ql- 


C 


(K,KA)+Q1+V1+W13] 


= (D02,DA02) 


Q2- 


C 


(K,KA)+Q2+V2+W24] = 


(D14,DA14) 


D02 d 


$ DA14 


DA02 d 


$ D14 


Ql- 


c 


(K,KA)+Q1+V1+W14] 


= (D03,DA03) 


Q2- 


c 


(K,KA)+Q2+V2+W23] = 


(D15,DA15) 


D03 d 


$ DA15 


DA03 d 


$ D15 


Ql- 


c 


(K,KA)+Q1+V1+W23] 


= (D04,DA04) 


Q2- 


c 


(K,KA)+Q2+V2+W14] = 


(D16,DA16) 


D04 d 


$ DA16 


DA04 d 


$ D16 


Ql- 


c 


(K,KA)+Q1+V1+W24] 


= (D05,DA05) 


Q2- 


c 


(K,KA)+Q2+V2+W13] = 


(D17,DA17) 


D05 d 


$ DA17 


DA05 d 


$ D17 


Ql- 


c 


(K,KA)+Q1+V1+W34] 


= (D06,DA06) 


Q2- 


c 


(K,KA)+Q2+V2+W12] = 


(D18,DA18) 


D06 d 


$ DA18 


DA06 d 


$ D18 


Ql- 


c 


(K,KA)+Q1+V2+W12] 


= (D07,DA07) 


Q2- 


c 


(K,KA)+Q2+V1+W34] = 


(D19,DA19) 


D07 d 


$ DA19 


DA07 d 


$ D19 


Ql- 


c 


(K,KA)+Q1+V2+W13] 


= (D08,DA08) 


Q2- 


c 


(K,KA)+Q2+V1+W24] = 


(D20,DA20) 


D08 d 


$ DA20 


DA08 d 


$ D20 


Ql- 


c 


(K,KA)+Q1+V2+W14] 


= (D09,DA09) 


Q2- 


c 


(K,KA)+Q2+V1+W23] = 


(D21,DA21) 


D09 d 


$ DA21 


DA09 d 


$ D21 


Ql- 


c 


(K,KA)+Q1+V2+W23] 


= (D10,DA10) 


Q2- 


c 


(K,KA)+Q2+V1+W14] = 


(D22,DA22) 


D10 d 


$ DA22 


DA10 d 


$ D22 


Ql- 


c 


(K,KA)+Q1+V2+W24] 


= (D11,DA11) 


Q2- 


c 


(K,KA)+Q2+V1+W13] = 


(D23,DA23) 


Dll d 


$ DA23 


DA11 d 


$ D23 


Ql- 


c 


(K,KA)+Q1+V2+W34] 


= (D12,DA12) 


Q2- 


c 


(K,KA)+Q2+V1+W12] = 


(D24,DA24) 


D12 d 


$ DA24 


DA12 d 


$ D24 



TABLE II: The table shows 48 possible combinations of the CMB subtracted WMAP Q band maps and 24 cross combinations 
to get rid of noise. 

Here we explain the steps followed to obtain the foreground power spectrum for Q band where the number of DA's 
are more than one. Similar steps are also followed for V and W band. 

• We smooth all the WMAP DA maps and the cleaned maps to one degree beam resolution. 

• Take the cleaned map CI from 4-channel combinations as given in table I and subtract it from Ql map. And 
similarly take the cleaned map C12 from 4-channel combinations and subtract it from Q2 map. We label the 
CMB subtracted maps as Di as given in table II. 
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• Cross correlate the DOl(Ql-Cl) with the D13(Q2-C12) as discuss in table II in equation (9) to get rid of noise 
term and obtain the foreground power spectrum. 

• Repeat the above steps for 24 possible combinations obtained by 4 channel cleaned maps given in table II and 
while cross correlate we choose pairs of Di where detectors are not common. Obtain the mean foreground power 
spectrum. 

Since, K and Ka band have only one DA, cross correlation isn't feasible. In these cases, estimate of contribution of 
noise bias is explicitly calculated from the number of observation and sensitivity of the detector and subtracted out. 
To avoid any kind of mixing, we subtract the 3 channel cleaned maps CI from let say K band and cross correlated 
with another 3 channel cleaned map C12 (where the detector used are not common) subtracted K map. So, the 
equation (9) in case of one DA will become, 

((a™ al (i) - a^ an (i))(af^ al *(j) - a% an *(j))) 

= {{Btaf m {i)af*{j))) + ((a^(i)a^(i))). 

= cf + cl f . (10) 



CMB subtracted WMAP K maps 


Cross combinations 


K- 


C(Q1V1W12) 


=D01 


K-C(Q2V2W34) 


=D13 


D01 (5 


$ D13 


K- 


C(Q1V1W13) 


=D02 


K-C(Q2V2W24) 


=D14 


D02 (5 


$ D14 


K- 


C(Q1V1W14) 


=D03 


K-C(Q2V2W23) 


=D15 


D03 (5 


$ D15 


K- 


C(Q1V1W23) 


=D04 


K-C(Q2V2W14) 


=D16 


D04 £ 


$ D16 


K- 


C(Q1V1W24) 


=D05 


K-C(Q2V2W13) 


=D17 


D05 $ 


$ D17 


K- 


C(Q1V1W34) 


=D06 


K-C(Q2V2W12) 


=D18 


D06 $ 


$ D18 


K- 


C(Q1V2W12) 


=D07 


K-C(Q2V1W34) 


=D19 


D07 (5 


$ D19 


K- 


C(Q1V2W13) 


=D08 


K-C(Q2V1W24) 


=D20 


D08 (5 


$ D20 


K- 


C(Q1V2W14) 


=D09 


K-C(Q2V1W23) 


=D21 


D09 (5 


$ D21 


K- 


C(Q1V2W23) 


=D10 


K-C(Q2V1W14) 


=D22 


D10 $ 


$ D22 


K- 


C(Q1V2W24) 


=D11 


K-C(Q2V1W13) 


=D23 


Dll £ 


$ D23 


K- 


C(Q1V2W34) 


=D12 


K-C(Q2V1W12) 


=D24 


D12 


$ D24 



TABLE III: The table shows 24 different combinations of the DA maps for K band and 12 cross combinations to get rid of 
noise. There is a corresponding set for the Ka band. 

The second term in right hand side Cf of equation (10) can easily be estimated using the relation, 

N g?(l°) 4tt ^ a\ 
1 Bf(K) N% ix ^ N obs - 

By substituting the value of back in equation (10) and subtracting it from the left hand side gives the 1° beam 
smooth foreground power spectrum (Cf). In general, foreground power spectra are expressed in terms of antenna 
temperature which can easily be converted from thermodynamic temperature using the table IV. Henceforth, all the 
results and plots of Cf is expressed in terms of antenna temperature. 

To get the foreground power spectrum outside the KQ85 mask, we mask the difference map Di's maps with combine 
KQ85 and WMAP5 point source mask and then smooth to 1 degree before cross-correlating the Di's where no common 
DA/detector are present to get rid of noise. Additional smoothing of 1° is applied outside the KQ85 mask to compare 
the foreground power spectrum with MEM maps. 

We can calculate the r.m.s temperature, AT rms , of the foregrounds using the relation, 

(AT rms ) 2 — "V^ — Cf \ 

where AT rrns is expressed in terms of antenna temperature. We obtain the frequency dependent AT rms for the 
five WMAP frequencies. We can model the rms foreground power in terms of three major galactic foregrounds - 
synchrotron, free-free and dust emission. We fit the AT rms at the five frequencies to obtain the full sky synchrotron 
spectral index, 



(11) 



6 



Frequency 
(in GHz) 


Conversion Factor 

gM 


Mean <jo 
(in mK) 


Mean FWHM 
(in degrees) 


23 


0.9867 


1.436 


0.807 


33 


0.9723 


1.470 


0.624 


41 


0.9581 


2.197 


0.4775 


61 


0.9095 


3.133 


0.326 


94 


0.8012 


6.538 


0.2038 



TABLE IV: Conversion factor from thermodynamic to antenna temperature where AT a — g(v)AT; g(v) — [x 2 e x /(e x — l) 2 ]; 
x=his/k B T ; T =2.725K. 




Multipole(£j Multipole 

FIG. 1: Top Left Panel: The model independent estimate of angular power spectrum for combined foregrounds at WMAP 
frequency starting from K band to W band. In this case Cf is smoothed by 1 degree beam and expressed in terms of antenna 
temperature. Top Right Panel: Plot of relative difference in estimated power (Q^em _ qMI^jqMem w ^h the multipole /, 
where 'MI' stands for the foreground power spectrum obtained from model independent analysis. Foreground power spectrum 
at K and Ka band is consistent with the WMAP Team. But for Q, V and W band we obtain slightly less power compared 
to MEM method. Bottom Left Panel: The angular power spectrum of combined foregrounds outside KQ85 mask at WMAP 
frequencies starting from K band to W band expressed in antenna temperature. Bottom Right Panel: Plot of relative power 
difference (cf /IEM — Cf^ 1 ) / Cf dEM with the multipole I clearly shows that at low multipole MEM underestimates the foreground 
power outside KQ85 mask. 

where (3d = 1.8, /3f = —2.14 are taken as a constant parameters. 

As a consistency check, we compare our results with that obtained by the WMAP team using Maximum Entropy 
method (MEM) analysis. We find that our methods detects marginally lesser foreground power spectrum for Q to 
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Frequency 
(in GHz) 


RMS Temperature 
using MI Analysis 
(in mK) 


RMS Temperature 
using MEM maps. 
(in mK) 


23 


2.795 


2.801 


33 


1.210 


1.213 


41 


0.731 


0.736 


61 


0.310 


0.314 


94 


0.198 


0.201 



TABLE V: Comparison of rms foreground power obtained from Model Independent analysis and MEM method. The rms 
foreground power matches closely with MEM method with a minor deficient seen in the Model Independent case. 



10.00 F 




100 

Frequency, v (in GHz) 

FIG. 2: The total foreground emission (black dots) spectra obtained using model independent analysis compared to sum of 
foreground components (deep blue line connecting black dots) over the full sky. The average synchrotron spectrum from K to 
W band is -2.6 assuming the free- free spectral index (3f = —2.14 and dust spectral index (3d = 1.8 are constant parameters. 
The average synchrotron spectral index from K-Ka band and Ka-Q band are -3.0 and -2.91 respectively. 

W band as compared to WMAP MEM method over the full sky. K and Ka band is quite consistent as shown in 
figure(l). The excess power in WMAP MEM method comes from the low multipoles. To estimate the foreground 
power outside KQ85 region, MEM maps are first mask with KQ85 mask and then smoothed to 1° beam resolution 
to get foreground power spectrum of effective smoothing of 1.414°. We found that MEM method underestimates 
the foreground power at low multipoles outside KQ85 mask. The advantage of model independent analysis is that 
foreground power spectrum is not resolution limited. For each band, we can estimate the foreground power spectrum 
up to multipoles limited by the beam of that band. 

A. Verifying the Model Independent Method using Monte-Carlo Simulations 

We carry out a set of Monte-Carlo simulations to estimate error on the rms foreground temperature obtained by 
the model independent analysis. A set of 1000 simulations of CMB maps are made using the HEALPix software for 
A 7 ' s ^ e =128 with proper beam-width and noise properties for each frequency channel of WMAP 5yr. We choose to 
downgrade WMAP maps to AT s ^ e =128 for our 1000 simulations because foregrounds are important at low multipoles. 
We use the publicly available Planck Sky Model (PSM) version 1.1, to generate the diffuse foreground maps. We 
estimate 48 (from 4 channel cleaning) and 24 (from 3 channel cleaning) cleaned maps from each realizations as 
discussed in section II. We subtract them from the input maps and then follow the steps as discussed in section III. 
Finally the rms temperature of foreground from 1000 realizations are recovered and compared with the input values. 
The mean and variance of recovered foreground rms values are given in table VI. For 1000 simulations, we get 72000 
cleaned maps for Nside=128. The simulations shows that our model independent foreground analysis recovers the 
foreground r.m.s power very well. 

While comparing the PSM recovered rms temperature of foregrounds from simulations with that obtained by MEM 
analysis and our analysis, we note that the PSM foreground templates significantly underestimate the level of galactic 
contamination. The PSM generated foreground power spectrum for K band is much lower than the MEM generated 
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MI K-band MEM K-barid Difference K-band 




FIG. 3: The first column is the combined foreground maps for K, Ka, Q, V and W maps from top to bottom obtained by model 
independent analysis. The second column represents the combined foregrounds maps for K, Ka, Q, V and W bands obtained 
by MEM method and the third column is MEM minus model independent foreground maps. Model independent foreground 
maps for K and Ka are obtain from 3 channel clean maps whereas the model independent foreground maps for Q, V and W 
band are obtained from 4 channel clean maps. 



power spectrum. But the differences decrease with the decrease of level of foreground contamination and for W band 
it matches with the MEM templates. This result has also been noted by other groups [12]. 
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Frequency 
(in GHz) 


ATrms (in fiK) 
Input PSM Template 


ATrms (in fiK) 
Extracted using Our Analysis 


23 


1440.73 


1440.64 ± 3.98 


33 


628.51 


627.67 ± 3.86 


41 


396.01 


396.29 ± 3.94 


61 


206.85 


207.23 ± 3.62 


94 


181.71 


182.27 ± 2.97 



TABLE VI: Comparison of input foreground rms temperature (expressed in mK) with the recovered foreground rms temperature 
obtained by Monte- Carlo simulations. 



IV. ESTIMATION OF SYNCHROTRON SPECTRAL INDEX OVER DIFFERENT REGIONS OF SKY. 

It is well established that spectral index of synchrotron emission varies significantly over the sky. We estimate 
the spectral index of synchrotron emission over different regions of the sky using the model independent foreground 
estimation. The regions are defined by the 192 coarser pixels of HEALPix Pixelization at Nside=4. Synchrotron 
spectral index can be calculated easily by knowing the frequency dependence of AT r ^ s . By definition AT F ms is given 

by, 

AT/ ms = |^f>rw} , (12) 



1 = 1 



where the index 'F : represents the combined foreground emission coming from dust emission, free-free emission, 
synchrotron emission and N p i x denotes the number of pixels at resolution Nside=4 in the given region of the map. 
AT r ^ s is generally expressed in terms of antenna temperature. 

Using equation (7), the cleaned map in real space can be written as a sum of CMB map and residual noise as, 

AT clean = + ^RN ^ 

Multiplying AT clean with 1° resolution beam B(n.h') and subtracting from equation (4), we get 

AT R (n) = J (AT F (n) + AT RN (n))B(n.h f ) dh' + AT N (h) 
= AT ,F (n) + AT N '(h), 

where, AT R (n) denotes the map of foreground plus residual noise. The beam smoothed foreground map can be 
defined as, 

AT fF (n) = j ' AT F (n)B(n.n')dn'. 

We can define a quantity AT r F ^ ns as, 

N pi x 

(AT r R ms f = —— J2 (AT?(n)AT b R (n)), (14) 

where the index 'a' and '6' represents the two independent detectors whose noise are uncorrelated. Using the above 
relation, (AT^J 2 can be written as, 

(AT r R ms f = J- g {AT'/ (n)ATf (n)> + (AT a N ' (h)AT b N ' (n)) 

= (AT r F ms f + (AT r N ms )\ (15) 



The assumption in the above calculation is that foreground and noise are uncorrelated. Since the noise for two 
independent detectors are uncorrelated, the second term vanishes and we retain only the foreground rms power. The 
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foreground rms temperature for the frequency channel with more than one DA's can be calculated as, 

N p i x 



AT r F ms = AT r R ms = \^-£ (&T?(n)AT b R (h)} \ . (16) 



For Q to W band, where the number of DA's is ranging from two to four, the foreground rms temperature can be 
calculated using the above equation. For example to calculate the quantity AT r ^ s for Q-band, the steps are as follows: 

• Smooth all the input WMAP DA maps and the cleaned maps to one degree beam resolution. 

• Subtract the cleaned map CI obtained by 4 channel combinations from Ql band and similarly subtract cleaned 
map C12 obtained by 4 channel combinations from Q2 band. The CMB subtracted maps are label as Di given 
in table II. 

• Take the cross product of DOl(Ql-Cl) with the D13(Q2-C12) over the region of sky defined by a single pixel 
of Nside=4 and take the sum over all the pixels excluding the pixels covered by WMAP5 point source mask as 
described in equation (14). 

• Repeat the above steps for 24 possible combinations produced by 4 channel cleaned maps and took the ensemble 
mean of it. As the noise properties of two difference maps Dl and D13 are uncorrelated, the second terms cancels 
out of equation (15). 

But for K and Ka band where only one DA is present, we subtract 3 channel cleaned maps CI from K band and 
subtract 3 channel cleaned map C12 from K band and took the cross product over the region of sky to get rid of 
residual noise coming from model independent analysis of CMB power spectrum. The equation (15) for the case of 
frequency channels with one DA becomes, 



where AT^ ns can be defined as, 



(AT rms ) (AT rms ) + (AT rms ) , 



(AT r N ms ) 2 = -}-Y^(AT a N (n)AT b N (n)). 



= 1 



which has been estimated using 1000 simulations of noise maps for K and Ka band smoothed to 1 degree beam 
resolution over the given region of the sky. The rms foreground temperature for K and Ka band can be calculated 
using the relation, 



AT 

rms 



N, 



— J^(AT R (n)AT b R (n))\ - AT r N ms . (17) 



= 1 



The AT r ^ s obtained is in terms of thermodynamic temperature which is converted to antenna temperature using the 
table IV. For each region, the frequency dependence of AT^ ms is fitted using the relation, 

where /3 d = 1.8, f3f = —2.14 are taken as a constant parameters. We fit AT r ^ s vs frequency to calculate the 
synchrotron spectral index at each region of the sky. To break the degeneracy, we use the Haslam map as a input 
template to increase the degree of freedom while computing the synchrotron spectral index. The rms noise level of 
Haslam map is 0.5K at high latitudes and 0.7K at low latitudes. Since the Haslam map isn't corrected for noise so 
the resulting spectral index variation puts an upper limit on the value of average (3 S . The spectral index variation at 
each pixel is shown in figure (4) as a sky map where we replace each HEALPix pixel with a colored circle which fits 
inside the pixel. We exclude the galactic region between ± 5° latitude in our analysis. 



11 




FIG. 4: Synchrotron spectral index variation over different positions of the sky obtained by model independent foreground 
analysis. The spectral index behaviour clearly shows that f3 s is -3.5 at high latitudes and -2.5 close to the galactic plane which 
is consistent with WMAP Team. For example the fit for 3 regions are shown in the plot. The region 126 (latitude range [10° < 
£ < 29.5°], longitude range[248° < b < 269°], synchrotron spectral index, \j3 a = -2.82]), region 137 (-59.1° < £ < -30.7°, 
23.1° < b < 44.5°, p 8 = -3.70), region 165 (-53.3° <£< -30.3°, 248° < b < 267°, f3 s = -3.33). 



V. DISCUSSIONS 

The estimation of foreground power spectrum from WMAP is carried out in a self contained method without 
using any extra information at any other frequencies other than WMAP frequencies. This work can be considered 
as an comprehensive approach that estimates both the CMB power spectrum and the foreground power spectrum 
simultaneously in a model independent approach. The method described in this paper is unbiased, we established 
through Monte-Carlo Simulations. Table VI shows that the recovered foreground rms power is very close to input 
foreground rms power we put using the PSM template. Interestingly we find that MEM method overestimates the 
foreground power close to galactic plane and underestimates the foreground power at high latitudes relative to our 
estimates. But for the full sky, the mean rms foreground power using model independent method is close to MEM 
method. We find that the average synchrotron spectral index from K to W band is f3 s = —2.6 over the full sky. The 
behaviour of synchrotron spectral index of (3 S = —3.5 at high latitudes and f3 s = —2.5 close to the galactic plane 
is consistent with WMAP Team. The advantage of this method is the cross-correlations takes care of the residual 
noise which remains after cleaning the map. For the upcoming PLANCK mission, this method of foreground power 
estimation will be even more promising since there are huge number of cross-combinations available due to large 
number of detectors of it and greater frequency coverage of PLANCK. 
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